Penetrating intrinsic charm: evidence in data 
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Nuclei are transparent for a heavy intrinsic charm (IC) component of the beam hadrons, what leads 
to an enhanced nuclear dependence of open charm production at large Feynman xf. Indeed, such 
an effect was observed recently in the SELEX experiment [lj. Our calculations reproduce well the 
data, providing strong support for the presence of IC in hadrons in amount less than 1%. Moreover, 
we performed an analysis of nuclear effects in J production and found a similar, although weaker 
effect. 

PACS numbers: 13.85.Ni, 14.40. Lb, 14.20.Lq, 14.40.Pq 
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All data available so far for particle production at for- 
ward rapidities in high-energy hadron-nucleus collisions 
exhibit a similar trend for enhanced nuclear suppression 
at larger Feynman xp (e.g. see the collection of data 
in 0]). However, the SELEX experiment has released 
recently data for charmed hadron production in hadron- 
nucleus collisions [lj] , which demonstrate an opposite be- 
havior: the nuclear ratio is rising, rather than falling, at 
large x f ■ Here we interpret this effect as a manifestation 
of the projectile intrinsic charm (IC) component, which 
does not attenuate propagating through the nucleus. So 
it is filtered out and relatively enhanced by the nucleus. 
Such a signal of IC can also be observed in data for nu- 
clear dependence of J/^> production. 

Charmed hadron production. First of all, one should 
understand why produced hadrons are suppressed by nu- 
clei, as one can see in the SELEX data [l[ for charmed 
hadrons, as well as in data for other hadronic species 
According to the Mueller-Kancheli theorem [l,[4j], inclu- 
sive cross sections of particle production, hA —> h'X, 
should have no shadowing. This conclusion is based on 
the Abramovsky-Gribov-Kancheli (AGK) cutting rules 
Q with an extra factor which is equal to the number of 
cut Pomerons. In the particular case of hadron-nucleus 
collisions one can expand the inclusive cross section over 
the number of collisions n co u, with the Glauber weights 
and with an extra factor, n co n, compared to the same 
expansion for the total cross section. This is because the 
final hadron h! can emerge from any of the multiple in- 
elastic collisions. Then, instead of the A 2 I Z characteristic 
behavior of the total cross section, one obtains a linear, 
oc A, rise of the inclusive cross section. 

This prediction has been indeed confirmed in hard re- 
actions having small cross sections, like the Drell-Yan 
process, production of hadrons and direct photons with 
high transverse momentum, etc. However, in soft pro- 
cesses one should expect the linear A dependence to be 
broken due to saturation of the unitarity bound. Namely, 
hadron production is enhanced by multiple interactions 
only if the cross section is far from the unitarity bound. 
Otherwise, the production processes on different nucle- 



ons shadow each other and the cross section remains es- 
sentially unchanged by multiple interactions, an effect 
known as Landau-Pomeranchuk principle (see in [||). 
Only the transverse momentum distribution is modified, 
the called color glass condensate [7| . It is also related to 
coherence in particle production from multiple sources: 
when the phase space of produced particles is densely 
packed up to some transverse momentum Q s , the inclu- 
sive cross section saturates with number of collisions, i.e., 
the inclusive cross section is independent of number of 
collisions. Then, one should suppress the extra factor 
n co u in the Glauber expansion, and eventually arrive at 
a shadowed inclusive cross section oc A 2 / 3 . 

Moreover, at large Feynman xf 1 the nuclear sup- 
pression is unavoidably getting stronger due to the dis- 
sipation of energy in multiple interactions in the nu- 
cleus, and the restrictions imposed by energy conserva- 
tion. Eventually, near the kinematic limit the suppres- 
sion reaches the maximal possible magnitude correspond- 
ing to the A 1 / 3 dependence of the cross section. One 
can interpret this mechanism as excitation of higher Fock 
components in the beam hadron by multiple interactions 
in the nucleus. The more constituents a Fock state has, 
the steeper is the fall of the parton distribution func- 
tion at x — > 1. Thus, each of the multiple collisions, but 
the first one, supplies an additional suppression factor 
S(xf) @i an d summing up the modified Glauber expan- 
sion over number of collisions, one gets 
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Here the effective atomic number A e ff{xp) — 
a A (x f) / cr N (x f) is the ratio of the inclusive cross 
sections. The nuclear thickness function 7a(&) = 
dz pA(b, z) is the integral of the nuclear density 
Pa(o, z) along the trajectory at impact parameter b. 

The suppression factor S(xf) should vary from 1 at 
xp = (see above) down to vanishing as 1 — xf at 
xp — > 1 0, HI- Here we rely on the simplest form 



S(xf) = 1 — Xf, which describes well a bulk of data 
for production of light hadrons @. The effective ab- 
sorption cross section a a b s , which controls the number 
of collisions, in the Glauber approximation would be just 
the inelastic hadron-nucleon cross section. However, the 
Gribov inelastic corrections Q make the nuclear medium 
more transparent and lead to a considerable reduction of 
a a bs compared to [l(| • Similar to Ref . Q , we evalu- 
ated the effective absorption cross section by comparing 
the Glauber exponential attenuation with one calculated 
in the dipole approach (see Eq. (58) in [H|). For this 
evaluation we use the mean nuclear thickness for each 
nucleus. For E~ hyperon and pion we use the results of 
the SELEX experiment [3, a^ t /a^ t = 0.9 and 0.635 for 
h = £~ and pion, respectively. This procedure was ap- 
plied to copper and carbon separately and as function of 
energy, in accordance with data in Ref. For example, 
for copper at pi a b — 300 GeV we obtained a^ bs = 20.4 mb, 
a abs = 13-1 mb and af b ~ = 18.8 mb . Notice that with 
this value of a v ahs Eq. fll]) describes well the xf depen- 
dence of nuclear suppression observed for many hadronic 
species produced in pA collisions [2J . 

In Fig.[T]our calculations for the exponent a are shown 
by dashed curves in comparison with the SELEX data. 
While at small x f our results agree with the data, Eq. (JTJ) 
underestimates the value of a measured at large xf ■ This 
anomaly might signal a missed IC contribution, since 
IC in light hadrons carries the main fraction of its mo- 
mentum [L3l - flq ] and becomes a new source of charmed 
hadrons at large xf- Since the mechanism of fragmenta- 
tion of a light projectile quark is nuclear suppressed at 
large x f , the fragmentation of a fast charm quark origi- 
nating from the IC becomes more important. 

A peculiar feature of heavy quarks is their weak at- 
tenuation in a medium. The mean momentum fraction 
of an intrinsic heavy component remains large even in 
higher Fock components, provided that the heavy quark 
mass squared is considerably larger than the square of the 
hadron mass. For the same reason the medium induced 
perturbative energy loss is much less for heavy than for 
light quarks [13] ■ Thus, multiple nuclear interactions, 
which excite higher Fock states in the projectile hadron, 
do not affect much the IC momentum distribution and 
do not suppress the originated from IC charmed hadrons 
at large x f ■ This is the key observation, which makes the 
IC a plausible mechanism responsible for the anomalous 
nuclear dependence observed in the SELEX experiment. 

To proceed further with the interplay of the two mech- 
anisms of charmed hadron production, notice that the 
limiting behavior of the cross section at xp — > 1 is inde- 
pendent of presence or absence of IC, since it is controlled 
by the triple Regge formalism, 
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where t is the 4-momentum transfer squared, and aR c (t) 
is a charmed Regge trajectory which assignment depends 



FIG. 1: (a) The exponent o(xf) characterizing the A- 
dependence of the cross section as function of xf. The dashed 
curves show o(if) for the conventional (no IC) mechanism 
of light quark fragmentation calculated with Eq. (fl]). Solid 
curves show the effect of inclusion of IC as is given by Eq. (|4| . 
The curves from bottom to top correspond to proton, E~ and 
pion beams, respectively, (b) Same as above, but for the com- 
bined statistics for all beams. 



on the quantum numbers of h and h c . Only the struc- 
ture the h-h c -R c vertex depends on the production mech- 
anisms under discussion, as is illustrated in Fig. but 
the Reggeon remains the same. 

The ^-integrated cross sections Eq. ^ at xf — > 1 has 
the form da/dxF oc (1 — xp ) n , which is the universal 
behavior for both mechanisms. We denote the relative IC 
contribution (Fig. [5^, right) by S, correspondingly the rest 
(Fig. [2^ left) is I — 8. While these two mechanisms have 
the same xf dependence at xf — ► 1, the IC contribution 
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FIG. 2: a: Regge vertex h-h c -R c for two mechanism of h c 
production. X c and X q denote quark assembles with or with- 
out open charm, respectively, b: The vertex p- J/»l/-Reggeon 
for two mechanisms of J/\£ f production, qq annihilation (left), 
and either the perturbative glue-glue fusion, or the IC con- 
tribution (right). Xs q and X2 q are the three or two quark 
partonic systems, respectively; X g is a pure gluonic system. 



declines from this behavior at smaller xp. In particular, 
the IC contribution is expected to have a maximum at 
xp w (zic)/2, where (zic) — 0.75 is the position of 
the peak in the IC distribution function [18] . Thus, we 
assume the xp -dependence of h c production to be shaped 
as 
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where m = (zjc) nj (2 - (z IC )) = 5 n. 

The exponent n controlling the asymptotic behavior of 
charmed hadron production at xp — ¥ 1 can be expressed 
via the intercepts of the relevant Regge trajectories. For 
instance, n(7r~_-> D°,D~) = 1 - 2a_D»(0) + A = 3.5; 
ji(tt- -> D+,D°) = 1 - 2arj.(0) - 2a p (0) + 2 + A = 
4.5. The term A « 2a'(\t\) = 0.5 is a result of the 
^-integration. The intercept and slope of the D* tra- 
jectory was evaluated in [lj[ at od»(0) = —1, and 
On, = 0.5 GeV~ 2 . So, for the reactions measured in 
nfl with the pion beam we chose = 4. Similarly we 
estimated n p — 5.5 and ny, = 6.5. Then, we are in a po- 
sition to add up the two mechanisms of charmed hadron 
production in the effective atomic number 



A eff (xp) = (l-5) A eff {x F ) + 6x F 1 A, 



(4) 



where A e ff(xp) for the conventional mechanism of 
light quark fragmentation was calculated earlier (dashed 
curves in Fig. [IJ . All the parameters have been evaluated 
above, except 8 which we treat as a fitting parameter. 
However, the experimental errors at large xp are too big 
for a statistical fit. One can get a fair description of the 
data within the range of 8 = 0.5 — 0.9. The solid curves 
depicted in Fig. Q] correspond to 8 = 0.8. 

As far as 8 is known, one can try to estimate the IC 
probability Pic- The total cross section of charm pro- 
duction from IC can be presented as Picnic, where 07c 
is a part of the total inelastic cross section in which the 
IC component is resolved and freed. Comparing this with 
the above estimates one arrives at 
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Here B(x, y) is the Beta Function; ao is the factor in 
the ^F-dependent cross section of inclusive production 



of charmed hadrons, da(hp — > h c X)/dxp = <7o(l — xp) n . 
We use the data from the LEBC-EHS experiment [2(| for 
D~ , D°-meson production in pp collisions at 400 GeV at 
CERN SPS, which were fitted in [H with n = 5.4±1.2, in 
good agreement with our evaluation n p — 5.5. Summing 
up these two channels and doubling the result to account 
for D production, we get ao = 170 Production of A c 
from hadronization of a c-quark is strongly suppressed at 
xp — > 1 and can be neglected. 

The magnitude of 07c dependents on a model for 
the structure of IC. The smallest cross section, i.e., the 
largest estimate Eq. ([5]), would correspond to an IC of 
perturbative origin with a cc separation of the order of 
X/m c . In this case 07c ~ lmb, and the range of possi- 
ble IC weights would be P IC ~ 0.0009 - 0.008. If the IC 
had a nonperturbative origin, and the cc were sitting in a 
nonperturbative potential, the typical cross section 07c 
would be like for J/\f r [2l|, ao ~ 5mb. Correspondingly, 
Pic ~ 0.0002 - 0.002. This estimate is much smaller 
than the previous evaluation in Ref. [22j based on the 
EMC measurements of charm distribution in the proton. 

Charmonium production. The penetrating IC com- 
ponent may also show up in the charmonium produc- 
tion as an enhanced A-dependence at large xp. How- 
ever, in this case IC cannot dominate the endpoint re- 
gion of xp — > 1 where the valence quark contribution 
takes over. Therefore an anomalous nuclear dependence, 
if any, should be a weaker effect. Fig. [2Jd illustrates 
the competing mechanisms of charmonium production 
in terms of the Regge approach. The first vertex corre- 
sponds to qq annihilation, and the related Regge inter- 
cept is a R3q (0) = a N (0) = -1/2. 

In the second vertex in Fig. [3t>, all three valence quark 
of the proton (the system X 3q ) get a low fractional mo- 
mentum, while the main fraction of the proton momen- 
tum is carried by the cc. It does not make any difference 
whether this pair was created perturbatively, g — > cc, or 
preexisted in the proton as IC. The i-channel exchanges 
and related xp dependences at xp — > 1 are identical, like 
in the case of open charm production (see above). The 
color-octet cc pair converts to J/'J' via gluonic exchange 
X g . Thus, instead of a color triplet-antitriplet system, 
in this case a color octet-octet dipole, {X 3q } s — {X g } 8 , 
is exchanged. In l/N c approximation this octet dipole 
can be presented as two triplet dipoles 2q — q and q — q. 
This corresponds to the p-N Regge cut with the intercept 
a£ 35 (0) = ajv(0) + a p (0)-l = -l. 

The endpoint xp dependence of J/^ production, 
dajdxp ex (1 — xp) n *, correlates with the mechanism: 
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1 - 2a7v(0) + A^ = 4; n| c = 3 - 2a N {0) - 2a p {0) + 

AjA = 5. Here = 2a' N (|t|)^ = 2 for {\t\)^ = 1.1 GeV 2 
[24|. Thus, r& < ri^ c , so the qq annihilation, rather than 
IC, dominates at xp — > 1. However, at medium high xp 
data favor the (1 — xp) 5 dependence [25[ indicating the 
dominance of the cc mechanism. The latter includes pro- 
duction of a cc either by projectile gluons, or from the IC. 
Both contributions have the same end-point behavior at 
Xp — > 1, but the IC part declines at smaller xp, since it 
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FIG. 3: The exponent characterizing the A Q -dependence of 
the cross section of J/ty production as function of xf- The 
dashed curves show the result of the fit to the E866 data [23l ] 
by Eq. (0 and extrapolation to larger xf- The solid curve 
corresponds to inclusion of IC. 



peaks at (xic) =0.75. To comply with this restrictions 
we choose the following form of the IC part of the J/^ 
production cross section, 



with fixed parameters (5 = 0.8 and e = 0.04 (fitted to the 
absolute cross section [13]). A q e fj(xp) was calculated in 
analogy to what we have done for open charm (without 
IC, and with added c-quark shadowing calculated in [26|). 
(xf) was parametrized as In^A^^/A) = a + bxp + 
The result of the fit with a = 1.11, b = -0.56, 



cx 



eff 
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0.27 is plotted by dashed curve in Fig. [3] 



Now we can add into Eq. ([7]) a third term presenting 
the IC contribution, which according to Eq. © has the 

3n 5c 

form 8 x F * A. The final result plotted by solid curve 
in Fig. [3] does not contradict the E866 data, although it 
is difficult to say with certainty that the data need it. 
Notice that the E789 experiment on J/ty production [24| 
also observed some disorder in A-dependence at large xf, 
which might be a result of the IC contribution. 

Summarizing, the anomalous nuclear dependence for 
charmed hadrons produced at large xf gets a natural 
interpretation in terms of the projectile intrinsic charm 
which does not attenuate in a nuclear medium. Calcu- 
lations performed here are in a good accord with data, 
provided that the probability of IC in the hadronic wave 
function ranges from 0.1% to 1%. A similar contribu- 
tion of IC to J/vp production also causes an enhanced 
A-dependence at large xp, which also agrees with data. 



do 



ic 



3n<V 



dxp 



oc x F * (1 - x F ) n * 
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The mechanism of nuclear suppression of J/4 1 , which 
is still under debate, is not the main focus of this paper. 
Therefore we simply fit the data by the gluon fusion and 
quark annihilation mechanisms at 0.3 < xp < 0.7, where 
the IC contribution is expected to be small, 



A J J^{xf) oc 



1 — xf 



Al ff (xF) + (l-6)A? ff ( XF ), (7) 
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